Time-dependent transport of electrons through a photon cavity 
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We use a non-Markovian master equation to describe the transport of Coulomb interacting elec- 
trons through an electromagnetic cavity with one quantized photon mode. The central system is 
a finite parabolic quantum wire that is coupled weakly to external parabolic quasi-one-dimensional 
leads at t = 0. With a stepwise introduction of complexity to the description of the system and 
a corresponding stepwise truncation of the ensuing many-body spaces we are able to describe the 
time-dependent transport of Coulomb-interacting electrons through a geometrically complex central 
system. We take into account the full electromagnetic interaction of electrons and cavity photons 
without resorting to the rotating wave approximation or reduction of the electron states to two lev- 
els. We observe that the number of initial cavity photons and their polarization can have important 
effects on the transport properties of the system. The quasiparticles formed in the central system 
have a lifetime limited by the coupling to the leads and radiation processes active on a much longer 
timescale. 

PACS numbers: 73.23.-b, 78.67.-n, 42.50.Pq, 73.21.Hb 



I. INTRODUCTION 

During the past decade there has been increasing in- 
terest in exploring time-dependent quantum transport 
through open mesoscopic systems in a strong system- 
lead coupling regime.^ ^ Utilizing the tunable dynamic 
response of transient time-dependent transport enables 
development of switchable mesoscale electronic devices, 
in which the interplay of the mesoscopic system with ex- 
ternal perturbations plays an important role.l^HI^ 

In the weak system-lead coupling regime, the wide- 
band and the Markovian approximation are usually em- 
ployed, by neglecting the energy dependence of the elec- 
tron tunneling rate, as well as memory effects in the sys- 
tem, respectively.'^SH^ It is assumed that the correlation 
time of the electrons in the leads is much shorter than 
the typical response time of the central system. However, 
transient transport properties, which are intrinsically 
linked to coherence and relaxation dynamics, cannot gen- 
erally be described in the Markovian limit. One has 
to take into account the energy-dependent spectral den- 
sity in the leads, and an accurate numerical method for 
such a nonequilibrium transient time-dependent trans- 
port is desirable. A non-Markovian density-matrix for- 
malism involving the energy dependent coupling elements 
should be considered based on the generalized master 
equation (GME).'^^^^How to appropriately describe the 
carrier quantum dynamics under nonequilibrium con- 
ditions in reali stic device geometries is a challenging 
problem.™mi 

More recently, manipulation of the electron-photon 
coupled quantum systems in an electromagnetic cavity 
has become one of the key issues to be implemented 
in quantum information processing devices. Utilizing 



the giant dipole moments of intersubband transitions in 
quantum wells^^ enables researchers to reach the ul- 
trastrong electron-photon coupling regime.'^^''^ In this 
regime, the dynamical electron-photon coupling mech- 
anism has to be explored beyond the wide-band and 
rotating- wave appr oximat ions .^^^^^^ Nevertheless, the dy- 
namical time-dependent transport of Coulomb interact- 
ing electrons in a specified geometry through an electro- 
magnetic cavity with quantized photon modes remains 
unexplored. 

In the present work, we explore the electronic tran- 
sient transport dynamics of an open quantum wire placed 
in a lin early polarized electromagntic field created in 
a cavity.'^Il^ The wire is contacted to two quasi-one- 
dimensional (QID) semi-infinite leads. A bias voltage 
is suddenly switched on between the leads, along the 
wire. The whole structure is considered to be placed 
in a parpendicular homogeneos external magnetic field. 
We use the Nakajima-Zwanzig (N-Z) formalism to project 
the time evolution of the system onto the Hilbert space of 
the central element (the short wire) by tak ing trace with 
respect to the operators in the leads .^^^^ The transient 
transport properties will be explored by calculating the 
time-dependent total mean number of electrons and pho- 
tons, and the time-dependent total charge current from 
the left (L) lead to the right (R) lead. 

The paper is organized as follows. In Sec. II, we shall 
describe our theoretical model including the electron sys- 
tem in an electromagnetic cavity by coupling a many- 
level electron system with photons using a full photon en- 
ergy spectrum of a single cavity mode in the Fock space. 
The N-Z fram ework is utilized to describe the system- lead 
coupling.l^^^ In section III we investigate the dynamical 
transient transport properties. Concluding remarks will 
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be presented in Sec. IV. 

II. MODEL 

We start by describing a Coulomb interacting electron 
system in a closed finite system in an electromagnetic 
cavity with quantized photon modes of two different po- 
larizations. Later we shall open the system by coupling 
it to external QID leads in order to allow for transport 
of electrons through it. In the closed system the num- 
ber of electrons is constant, but in order to accomplish 
the opening up of the system we need to calculate its 
equilibrium properties for a varible number of electrons. 
In the equilibrium and transport calculation to follow 
A^SES will represent the number of Single-Electron States 
(SESs) used to build the Many-Electron States (MESs) 
of the system.^^^i^ 

A. Electron system in an electromagnetic cavity 

The central element, a short wire placed in the x-y- 
plane is in an external homogeneous classical magnetic 
field constant in time, B = 5z, introduced by the vector 
potential Aext, and cavity fields represented by a quan- 
tized vector field A. The many-body Hamiltonian of this 
closed system is 

+ i?coui + -f^EM ~~ J dry A - ^^^^^ J 

= y I ^ + \m*^lly'^ + iJcoui + -ffEM 

with i^coui the Hamiltonian for the Coulomb interaction 
of the electrons and i^EM the Hamiltonian of the cavity 
photons to be introduced later and the charge current 
density and charge density 

j = -^{V^^(^V^) + (^>^)V^}, P=-eV^V, (2) 
where 

7r= (p+^Aext). (3) 

The static magnetic field B = V x Aext together with the 
parabolic confinement of the q-lDES introduce a charac- 
teristic length = ^/h/{rrlH}^, with 17^ = uj'^-\-Qq and 
the cyclotron frequency oOc = eB/{m*c). The frequency 



characterizes the strength of the electron confinement 
in the ^/-direction. The finite parabolic quantum wire has 
length Lx and hard walls at x = ±Laj/2. 

In terms of creation and annihilation operators the 
Hamiltonian of the closed system takes the form 

Ho = Hem + ^Coul + ^e-EM 

= ^Eidldi + hwa^a + ^ ^{ij\Vcoui\rs)dldldsdr 

i ijrs 

+ £cJ24dj 9ij{a + a^} (4) 

+ (^) E 4d. { {a^a +l)+l{aa + a^a^) } , 

where the single-electron states (SESs) of the closed sys- 
tem are labeled with Latin indices, r, 5}, a is the 
destruction operator of one quantum of the single-mode 
cavity field with frequency cj, and di is an annihilation 
operator of the non-interacting single-electron state \i) 
with energy Ei. 

The electromagnetic cavity is a rectangular box 
{x,y,z) e {[-ac/2,ac/2] x [-ac/2,ac/2] x [-4/2,4/2]} 
with the finite quantum wire centered in the z = plane. 
The polarization of the electric field can be chosen along 
the transport direction, x, or perpendicular to it by se- 
lecting the TEoii or the TEioi mode, respectively. In 
the Coulomb gauge the vector field is then 

with the upper equation representing the TEqh mode 
and the lower one for TEioi- 

The effective dimensionless coupling tensor of the elec- 
trons to the cavity mode due to the linear term in A in 
Eq. ^ is 

= S /^r[^*(r){(e-7r)V^,(r)} (6) 
+ {(e-7r)V^,(r)}>,(r)], 

since we have introduced the characteristic energy scale 
Ec = eA^w(iw/c = for the electron-cavity photon 
coupling. In the calculation of the energy spectrum of 
the electron-photon Hamiltonian Q we will retain all 
resonant and antiresonant terms in the photon creation 
and annihilation operators and not use the rotating wave 
approximation, but in the calculations of the electron- 
photon coupling tensor ([6| we assume a^^Lx « and 
approximate cos(7r{x, y}/ac) 1 in Eq. ([5| for the cavity 
vector field A. 

After the construction of the Fock space with the 
SESs that will later be deemed as relevant to the elec- 
tron transport an exact numerical diagonalization is used 
to obtain the Coulomb interacting MESs {|/i)} with the 
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energy spectrum and the unitary transformatiorp2125J 



Im) = ^V^a\a). 



(7) 



The total Hamiltonian of the electrons and the cavity 
fields can then be written in terms of the interacting 
MESs 

Ho = ^ \li)E^{li\ + hhja^a 



I ^a'^a + + i (aa + a'^a'^) I (8) 



and its energy spectrum has to be sought in a Fock space 
constructed from the space of the Coulomb interacting 
MESs {|a)} and the Fock space of photons {|A/ph)} 



-EM- 



(9) 



The diagonalization of the electron-photon Hamiltonian 
Q yields 



<^/e-EM- 



(10) 



The general scheme here is to start the numerical cal- 
culations with a fairly large number of SESs, A/'ses^ and 
retain only a certain number of MESs in the energy range 
that is relevant to the transport later, because the total 
number of MESs, A^mes, may be too large for the sub- 
sequent calculations. We will call a further twist on this 
procedure: "A stepwise introduction of complexity to the 
model and a stepwise truncation of its many-body space" . 
First, we introduce the Coulomb interaction between the 
electrons, second, we truncate the huge many-electron 
space. Then we add the Hamiltonian of photons and the 
electron-photon interaction, and again undertake a trun- 
cation of the many-body space in order to have a number 
of Many Body States (MBSs) for which the transport cal- 
culation can be performed. The MBSs are eigenstates of 
the interacting electron-photon system. We stress here 
that the unitary transformations necessary between the 
different many-body spaces, Eqs ([7|) and (10), have to be 
completed before truncation. In actual numbers, in the 
transport calculations here for = £c ^ 0.1 meV we 
used 10 SESs resulting in 1024 MBSs. The Fock-space 
of Coulomb interacting electrons was then truncated to 
the lowest 64 states. The inclusion of the single photon 
mode was accomplished with 27 photon states and the 
ensuing 1728 photon-electron many-body states will be 
truncated to the 64 states lowest in energy before the on- 
set of the transport calculation. In order to present the 
energy spectra for the closed electron-photon system in 



Fig. [T] we use 200 electron states and 20 photon states 
for the much larger values of the coupling g^^. For the 
closed system we calculate the energy-spectra for each 
number of electrons present separately. We present the 
calculations for the closed system in a separate publica- 
tion with more details on the convergence of the calcula- 
tions and properties of the system for high values of the 
coupling constant. There we also discuss when the sec- 
ond order term in A in the electron-photon interaction 
is necessary for our model of a finite quantum wire. 
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FIG. 1. (Color online) The many-body energy spectra for 
Coulomb interacting electrons coupled to quantized cavity 
photon modes with the electric component polarized along 
the finite quantum wire (x-polarization, left panel), and per- 
pendicular to the wire (y-polarization, right panel) versus the 
electron-photon coupling strength — Sc- In the energy 
range shown there are states with no electrons (green), one 
electron (red), and two electrons (blue). B = 0.1 T, hQo = 1.0 
meV, hw 0.4 meV, 300 nm, m* 0.067me, k 12.4, 

the dielectric constant of GaAs. 

The many-body energy spectrum of the electron- 
photon states are shown in Fig. [l] for the two polar- 
izations, along the finite quantum wire (x-polarization) 
and perpendicular to it (y-polarization). The horizontal 
states (green) in Fig. [l] are states only with photons and 
no electrons that are thus independent of the coupling of 
the electrons and the photons. 

In order to gain further insight into the character of 
the many-body states in Fig. [l] we plot in Fig. [2] the 
number of electrons in each state |/i), a conserved quan- 
tity, and thus an integer. In the same plot we add the 
diagonal of the photon number operator A^ph = a^a in 
the appropriate space {|/i)}. The photon number is not 
a conserved quantity in the closed system, but the diago- 
nal of the operator should give us some indication of the 
strong mixing of the photon-electron states by the cavity 
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FIG. 2. (Color online) The number of electrons in the many- 
body state (red), and the diagonal element of the photon- 
number operator, N^]l (blue). B = 0.1 T, g^^ = 0.1 meV, 
hQo 1.0 meV, hw 0.4 meV, 300 nm, m* 0.067me, 

= 12.4. 



coupling. Indeed, we see a considerable photon content 
in many states, not only the few ones that contain no 
electrons (i.e. /i = 1, 2, 5, 9, 14, 21, 31, 44, 58). 

The spectra versus the coupling constant in Fig. 
[ijhave resemblance with the energy spectrum of a single 
electron versus the magnetic field 5 in a non-circular 
quantum dot !^^*^^ * The reason comes clear when the 
single-electron Hamiltonian for the electron is written in 
terms of lowering and raising operators and compared to 
the many-body e-EM Hamiltonian Q. As long as the 
A^-term is included the magnetic field in the dot, 5, and 
the couling to the photons, in our model play a sim- 
ilar role. The x?/-symmetry breaking is not only caused 
by the geometry of the central system here, but also by 
the polarization of the photon field. 

Here, we have coupled a many- level electron system 
with photons using the full electromagnetic coupling ([8| , 
and it is interesting to compare the results with spectra 
obtained for the Jaynes-Cummings model without the ro- 
tating wave approximation derived by Feranchuk et 
and Li et al.'^^ This comparison will be detailed in an- 
other publication concentrated on the properties of the 
closed system. 



B. System connected to leads 

The closed system of Coulomb interacting electrons in- 
teracting with the single electromagnetic cavity mode is 
coupled at t = to external leads acting as electron 
reservoirs. We use a formalism proposed by Nakajima 
and Zwanzig to project the time evolution of the system 
onto the central system by partial trac ing op erations with 
respect to the operators of the leads .'^^^ The coupling 
Hamiltonian is of the form 



i,l 



X'{t) j dq {T^ct,d, + (T^J*4c,j}, (11) 



where / = {L, R} referring to the left and the right lead, 
and {t) is the time-dependent switching function of the 
coupling. The semi-infinite leads are in the same perpen- 
dicular constant external magnetic field B as the central 
system. Their energy spectra are continuous bands and 
the integral in Eq. (11) represents a summation over the 
band index and an integral over the continuous "momen- 
tum" q from the band bottom to an appropriate band 
cut-off. The coupling tensor • of a single-electron states 
\q) in the lead / to states \i) in the system is modeled as 
a non-local overlap integral of the corresponding wave 
functions in the contact regions of the system, 1^^, and 
the lead l^^^o. 



rjiL 

The function 



drdr'[4{r')]*^f{r)gl{ry). (12) 



glg{r,r') = g^exp [-S[{x - x'f - Sl^iy - y') 



/\21 



X exp 



\E^-e\q)\ \ 



(13) 



with r G and G Q^i defines the 'nonlocal overlap' 
and their affinity in energy. 

The Liouville-von Neumann equation describing the 
time-evolution of the total system, the finite quantum 
wire, the cavity photons, and the leads 

ihW(t) = [H(t), W(t)l W{t < 0) = pLPRPs, (14) 

where W is the statistical operator of the total system 
and the equilibrium density operator of the disconnected 
lead / G {L, R\ with chemical potential jii is 



Pi 



Tr^{e-^(^^-/^'^^)}' 



(15) 



Commonly, the spectral density Jo(c^) is used to de- 
scribe the coupling of states in an open system to reser- 
voirs or leads. ^^^^ We do not use it explicitly here, but 
the spectral density for lead / with respect to the SES a 
in the central system 



5{E-e\q)) 



(16) 



is a convenient tool to demonstrate graphically the phe- 
nomenological coupling selected here, in Equations ([l2| 
and ([l3|. In order to do this we show first in Fig. [3] the 
single-electron energy spectrum in the leads and in Fig. [4] 
the spectral density J\ {E) and the probability density for 
each SES used to build the MESs of the system. We see 
in Fig. [4] from jumps in the spectral density that SESs be- 
longing to different subbands in the central system have 
different coupling strengths to the subbands in the leads 
as a result of their symmetry. The simple geometry of the 
central system here does not result in states with widely 
different localization character. 
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FIG. 3. The single electron energy spectrum of the leads. 
B = 0.1 T, hQl = 1.0 meV. 



The Liouville-von Neumann equation is projected on 
the central system of electrons and photons by partial 
tracing operations with respect to the operators of the 
leads. Defining the reduced density operator (RDO) of 
the central system 



ps{t) = TrLTrRl^(t), ps{0) = ps, 



(17) 



we obtain an integro-differential equation for the RDO, 
the generalized master equation (GME) 



psit) = -^[JTs.PsW] 



Tr, 



{ 



HTit'), f dt' [U{t - t')H^{t')U+{t - t') 
Jo 

, Uo{t-t')ps{t')U+{t-t')]]}, (18) 



with the time evolution operator for the isolated sys- 
tems of Coulomb interacting electrons interacting with 
photons, and noninteracting electrons in the leads, 
U{t) = exp{-i{He^Hcoui^HEM^HL^HR)t/h}, 
without the coupling HT{t). The time evolution of the 
closed isolated system of Coulomb interacting electrons 
interacting with the photons is governed by Uo{t) = 
exp {—iHot/h}. 



The equation for the RDO (18) can be cast into a cou- 
pled set of integro-differential equations for the matrix 
elements of the RDO by projecting it on the {|/i)}-basis 
in a similar manner as was done for the Coulomb inter- 
acting electrons in our former work without photonsi ^^*^^ ' 
The kernel of the equation of motion for the RDO (18) 



has been approximated to the second order for the lead- 
system coupling i^T(^), but due to the integral structure 
of the equation the coupling is present in the solution to 
higher order. 

It is convenient to express many operators of the sys- 
tem in a particular basis, i. e. like the basis of 
independent electrons, or in the Coulomb interacting 
basis before coupling to the photons {|/i)}. The time- 
dependent mean values of these operators are best cal- 
culated by employing a unitary transformation between 
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FIG. 4. (Color online) The spectral density fi{E) (left) 
and the probability amplitude for the corresponding single- 
electron eigenstate (right). B = O.IT, hQo = 1.0 meV, 
= 300 nm, = 0.25 meV, g^al^^ = 53.14 meV, 

S[^^al = 0.4916. 



the bases. The photon number operator A^ph = a^a is 
best constructed in the {|/i)e-EM}-basis, and the unitary 
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transform introduced in Eq. (|To| to the {|/ti)}-basis yields 



(19) 



:0, g=0.01 meV 



= Trs{/9(t)W+iVphW}. 



Just as a reminder, the unitary transform is performed 
before the basis is truncated to the size of the matrix 
representing the RDO in the {|/i)}-basis dictated by the 
computational effort in solving the GME (18). 



The time-dependent average current into the sys- 
tem from the left lead, and from the system into the 
right lead are calculated from the GME (p!8|) along direc- 
tions introduced in earlier publications. 



III. TRANSPORT PROPERTIES 

For a smooth coupling of the left and the right leads 
to the central system of electrons and photons we use a 
switching function (12) 



1 



I e {L,R} 



(20) 



with 



0.2 ps ^. We fix the temperature of the 



reservoirs at T = 0.5 K, and the coupling strength as 

7 3/2 

go^^w = 53.14 meV. We select a bias over the central 
system by specifying jul = 2.0 meV, and jj^r = 1.4 meV. 

The time-evolution of the total mean number of elec- 
trons (A^e(^)) is displayed in Fig. [s] for two values of the 
electron-photon coupling g^^. Initially, no electron is 
present in the central system, but an integer number of 
photons is specified with a given polarization. For the 
weaker coupling, = 0.01 meV, the charging of the 
central system seems to be rather independent of the ini- 
tial number of cavity photons present and their polar- 
ization. Well after the charging, {t ^ 300 ps), a slight 
difference in the total charge as a function of the ini- 
tial condition can be observed, but it is very similar to 
the difference seen for higher coupling, g^^ =0.1 meV. 
More interestingly, the charging is influenced by the ini- 
tial photon number and their polarization for the higher 
electron-photon coupling. The charging is slowed down 
by a higher number of x-polarized photons. The total 
number of electrons in the system does not exceed one as 
could be expected by comparing the values of the chemi- 
cal potentials in the leads and the energy spectral in Fig. 
[T] we are in the Coulomb-blocking regime. 

For the parameters chosen here describing the geome- 
try and the potentials describing the central system the 
aspect ratio of the system for the low energy range is 
such that states that can be associated with motion in 
the x-direction are more numerous than states describing 
motion in the ^/-direction as can be verified by a glance 
at the lowest 10 eigenstates shown in Fig. [4j We should 
also have in mind that the photon energy hu = OA meV 
< hfto = 1.0 meV, the characteristic energy for parabolic 
confinement in the ^-direction. 
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FIG. 5. (Color online) The total mean number of electrons 
{Ne{t)) as a function of time for — 0.01 meV (upper 

panel) and g^^ = 0.10 meV (lower panel). B = 0.1 T, hQo = 
1.0 meV, La: 300 nm, huj 0.4 meV, /xl = 2.0 meV, /j^r 
1.4 meV, = 0.25 meV, gl^al^^ = 53.14 meV, 6[^2al = 
0.4916, m* = 0.067me, and k = 12.4. 



A clearer picture of the charging phenomena can be 
obtained by observing not the total amount of charge, 
but the time-dependent charging of the individual MBSs 
that are included in our calculation. This information is 
presented in Fig. [gI for the weaker electron-photon cou- 
pling, and in Fig. [T] for the stronger one. 

We immediately notice that in the case of the weaker 
coupling. Fig. [6j the occupation of the MBSs is almost 
independent of the polarization of the cavity photons, 
but as could be expected a higher number of them ini- 
tially present promotes the occupation of higher energy 
states. The occupation of the MBSs for 4 initial photons 
is indicative of a system strongly out of equilibrium with 
some lower lying MBSs almost empty. 

In the case of the stronger electron-photon coupling the 
occupation of the MBSs is still similar to the results for 
the weaker coupling for the ^/-polarization, but a drastic 
change is visible for the x-polarization with more MBSs 
partially occupied. 

In order to analyze the effects of different number 
of photons on the charging of the central system we 
show first the time-evolution of the total mean num- 
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FIG. 6. (Color online) The mean number of electrons {Ne{t)) 
in a MBS for x-polarization (left) and ^/-polarization 
(right) as a function of time. Initially, at t = 0, there is 
no photon in the cavity (top panels), 1 photon (middle pan- 
els), or 4 photons (bottom panels). The initial number of 
electrons is zero in all cases. B = 0.1 T, ^™ = 0.01 meV, 
huj 0.4 meV, /xl 2.0 meV, /j^r 1.4 meV, hQo 1.0 
meV, = 0.25 meV, gl^al^^ = 53.14 meV, 6[^2al = 0.4916, 
La: 300 nm, m* 0.067me, and k 12 A. 



FIG. 7. (Color online) The mean number of electrons {Ne{t)) 
in a MBS \fl) for x-polarization (left) and ^-polarization 
(right) as a function of time. Initially, at t = 0, there is 
no photon in the cavity (top panels), 1 photon (middle pan- 
els), or 4 photons (bottom panels). The initial number of 
electrons is zero in all cases. B = 0.1 T, ^™ = 0.10 meV, 
huj = 0.4 meV, /xl = 2.0 meV, /j^r = lA meV, hQo = 1.0 
meV, = 0.25 meV, g^al^^ = 53.14 meV, ^1,2^1 = 0.4916, 
La: 300 nm, m* 0.067me, and k 12.4. 



ber of photons (A/'ph) for the two different values of the 
electron-photon coupling strength in Fig. [Sj In both 
cases the mean number of photons does not vary much 
for the ^/-polarization. Larger deviations are seen for 
the x-polarization with shorter period oscillations for the 
stronger coupling. 

Again, a better insight into the photon-electron dy- 
namics can be gained by observing the evolution of the 
mean photon number of each MBS which is presented 
in Fig. [9] for the weaker coupling and in Fig. [lO] for the 
stronger one. For both cases we notice the rapid decay of 
the photon number from their initial state if the number 
was initially 1 or larger. With reference and comparison 
to Fig. [2] we observe a fast build up of electron-photon 
many-body states with a considerable photon component 
and higher energy. The total number of photons does 
not change fast in the system, but the introduction of 
electrons to the system through the leads results in a 
fast redistribution of the photons into many-body states 
of quasiparticles, or in other words the preexisting ini- 
tial photons participate in and facilitate the build up of 
many-body states with higher energy. The slow change 
in the total photon number in Fig. |8] reflects phenom- 
ena of radiation absorption and emission while the rapid 
change with time in the occupation in Figs [6] and [7| indi- 
cates the creation of quasiparticles with definite electron 



and photon components. We thus observe phenomena on 
at least two different time scales. 

We notice again (see Fig. [9| that for the weak electron- 
photon coupling there is not a large difference between 
the two polarizations except in the case of no initial pho- 
ton. 

It is well known that the GME in the approximation 
used here can lead to nonphysical negative probability 
for occupation of individual many-body levels if the cou- 
pling to the reservoirs is too strong. We have avoided 
this problem by keeping the coupling low enough. It is in- 
teresting to note here that it is our experience that after 
coupling to the photon system such a negative proba- 
bility will first turn up in the probability for a photon 
occupation before it would be seen in the probability for 
the occpation of an electron. 

The total current from the left lead into the system 
and from the system to the right lead is displayed in Figs 
11 and 12 for both polarizations. The irregular short 



period oscillation seen in the current stems from addi- 
tion of partial currents per each many-body state with 
simple almost harmonic oscillations of various short pe- 
riods characteristic for each component. For the weak 
coupling, = 0.01 meV, and no photon present at 

t = the left- and right currents are almost identical 
and very close to the current for the electron system with 
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FIG. 8. (Color online) The total mean number of photons 
(Nph) as a function of time for = 0.01 meV (upper panel) 
and — 0.10 meV (lower panel), x-polarization (red), and 
^/-polarization (blue). B = 0.1 T, /iQo = 1.0 meV, huj = 0.4 
meV, /iL 2.0 meV, ^lr 1.4 meV, 300 nm, 0.25 



\T I 3/2 

meV, g^a^ -- 
and K = 12.4. 



53.14 meV, Si = 0.4916, m* = 0.067me 



no coupling to cavity photons. This changes slightly for 
some photons initially present in the system, especially 
for t ^ 300 ps when the system is approaching a steady 
state. Here, we should mention that the photon energy, 
hbj = 0.4 meV, was selected to be smaller than the char- 
acteristic confinement energy in the ^/-direction and not 
to be in resonance with any Coulomb interacting MESs. 

In the case of a stronger electron-photon coupling, 
^EM = 0.1 meV, the charging of the central system is 
attenuated by the presence of x-polarized photons at 
t = 0. We are describing a central system here in the 
absence of any potential that can give its eigenstates dif- 
ferent localization character, see Fig. |4j we do thus not 
expect any simple phenomena of photo enhanced con- 
duction. The photon energy hw = 0.4 meV is far from 
being in resonance with the characteristic energy for the 
7/- confinement, but it is much closer to the characteris- 
tic energy for the x-direction. Inspection of the charge 
distribution of, or the occupation of many-body states 
l/i) in Fig. [rjand the distribution of the photons in the 
same states in Fig. 10 demonstrates a strong correla- 
tion between the states and that the initial photons have 




lM-)o = |2) = lOeMph). g =0.01 meV, x-pol 



lM.)o = |2) = lOeMph). g^ =0.01 meV, y-pol 
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FIG. 9. (Color online) The mean number of photons (A^ph(t)) 
in a MBS for x-polarization (left) and ^-polarization 
(right) as a function of time. Initially, at t = 0, there are no 
photons in the cavity (top panels), 1 photon (middle panels), 
or 4 photons (bottom panels). B = 0.1 T, g^^ = 0.01 meV, 
huj 0.4 meV, /xl = 2.0 meV, fiR = lA meV, hQo = 1.0 
meV, = 0.25 meV, gl^al^^ = 53.14 meV, 6[^2al = 0.4916, 
Lx — 300 nm, m* — 0.067me, and k — 12.4. 



caused the electrons to be distributed into many states. 
In a linear response situation which we do not have here 
we would say that the photon scattering of the electrons 
reduced the conduction. Here we have to say that the 
stronger photon-electron interaction for the case of the 
x-polarization attenuates the charging of the central sys- 
tem. Observation of the partial currents through the 
many-body state, not shown here, tells the same story. 
There are more states contributing to the charging, but 
they all bring less charge into the system. To understand 
this situation it is good to have in mind the similarity 
between the cavity photons here and confined phonons. 



IV. SUMMARY 

In this article we have taken the first steps to describe 
the transport of Coulomb interacting electrons through 
a photon cavity taking into account the geometry of the 
central system and the leads and allowing for a strong 
coupling to the photon mode. The central system is a 
parabolic quantum wire of finite length. Its aspect ratio 
and confinement characteristics make excitations for low 
values of the external magnetic field easier in the trans- 
port direction, the x-direction, than perpendicular to it, 
in the ^/-direction. By selecting the energy of the photon 
mode below the characteristic confinement frequency in 



9 



Mo = I1) = lOelOph). g™=0.1 meV, x-pol |n)o = |1) = lO^ flp^), g™=0.1 meV, y-pol 




FIG. 10. (Color online) The mean number of photons (A^ph(t)) 
in a MBS for x-polarization (left) and ^/-polarization 
(right) as a function of time. Initially, at t = 0, there are no 
photons in the cavity (top panels), 1 photon (middle panels), 
or 4 photons (bottom panels). B = 0.1 T, g^^ = 0.10 meV, 
hw = 0.4 meV, fiL = 2.0 meV, /j^r = 1.4 meV, hQo = 1.0 
meV, = 0.25 meV, g^al^^ = 53.14 meV, 6[^2al = 0.4916, 
Lx — 300 nm, m* — 0.067me, and n — 12 A. 

the ^/-direction we demonstrate that the transport be- 
comes very dependent on the polarization of the cavity 
mode and the number of photons initially present in the 
system at t = 0, before the electrons enter it. Gener- 
ally, a higher initial photon number reduces the transient 
charging of the system with electrons. This effect is en- 
hanced by increased coupling between the electrons and 
the photons. The largest reduction is found for several 
x-polarized photons in the system. The photon energy is 
then high enough to disperse the electrons into numerous 
excited states. The system does only contain delocalized 
electron states and this behavior has resemblance to re- 
duced conductance due to phonon scattering of electrons. 
It has though to be kept in mind that we are observing 
a transient behavior here in the charging phase of the 
central system. 

We observe dynamics on two different time scales here, 
the fast dispersion of electrons and photons into excited 
states initially seen in the transient regime, especially 
when analyzed for each many-body state, see Figs |6j [7| 
[9j and [lOj and the slow decay or gain caused by radi- 
ation phenomena and reflected by the total number of 
photons in Fig. [Sj Gonnected to the issue of the two 
time scales it is interesting to notice that in the case of 
no photon present in the system initially only many-body 
states with low energy and almost vanishing photon con- 
tent are occupied by the electron entering the system. In 



Nph=0, g'='^=0.01 meV 
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t(ps) 

Nph=4, g^'^=0.01 meV 
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t(ps) 

FIG. 11. (Golor online) The total current from the left lead 
(L), and into the right lead (R) as a function of time for 
^EM _ Q Q-j^ meV. Initially, at t = there are no photons in 
the cavity (top panel), or 4 photons (bottom panel). B = 0.1 
T, HVLq 1.0 meV, 300 nm, huj 0.4 meV, iil 2.0 

meV, iiR 1.4 meV, 0.25 meV, giaj'^ 53.14 meV, 

^1,2^1 0.4916, m* 0.067me, and i<i 12.4. 

this situation the total photon content of the system only 
changes on the radiation time scale. 

In the leads we have only electrons and in the central 
system quasiparticles with photon and electron content. 
Their lifetime is determined by the coupling to the leads 
and the time scale of radiation processes. The original 
coupling Hamiltonian of the leads and the central sys- 
tem, i^T, describes the entry or exit of noninteracting 
electrons to the central system from the leads. The two 
unitary transformations of i^x, fii"st to the basis of inter- 
acting MESs and the second one to the electron-photon 
many-body basis guarantee that the coupling between 
the leads and the central system describes how electrons 
leave or enter the leads, and conversely how the number 
of quasiparticles or photons changes in the system. This 
stepwise introduction of complexity to the model, concur- 
rent stepwise truncation of the ensuing Fock-space, and 
unitary transformation of the Hamiltonian at each step 
to the appropriate basis enables us to attack a problem 
that otherwise requires too large Fock space for numerical 
calculations. This proceedure can still be streamlined in 
order to describe systems with stronger electron-photon 
coupling or more complex geometry. 
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or "stiffer" in the y- than the x-direction. We thus do 
not see much change in the properties of the system (see 
Fig. 



11 



and 
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12 ) varying the initial photon number or the 
couphng strength to the photons for the parameter range 
explored here. For the low coupling strength, = 0.01 
meV there is not much difference either between the sys- 
tem properties for the two polarizations, but it is interest- 
ing to notice how similar the charging per MBS for the 
stronger coupling, = 0.1 meV, and ^/-polarization 

(see the right panels of Fig. [t]) is to the charging at the 
weaker coupling and x-polarization (see the left panels 
of Fig. |6|. The same analogy is found for the photon 
component in the MBS (compare the right panels of Fig. 
10 to the left panels of Fig. |9|. This analogy hints at 
a highly anisotropic "effective coupling" of the electrons 
and cavity photons for the geometry selectrd here. 

The strong coupling of the electrons to the photon field 
leads to polarization effects in our system that require a 
large functional basis for the electron states of our sys- 
tem to describe correctly. This together with other im- 
portant properties of the closed system will be discussed 
elsewhere^ and is the limiting factor here in how strong 
coupling we can describe in our present system and still 
manage observation of the time evolution according to 
the GME. 



FIG. 12. (Color online) The total current from the left lead 
(L), and into the right lead (R) as a function of time for 
^EM _ Q meV. Initially, at t = there are no photons in 
the cavity (top panel), or 4 photons (bottom panel). B = 0.1 
T, hQo 1.0 meV, 300 nm, hu 0.4 meV, fiL 2.0 

meV, fiR 1.4 meV, 0.25 meV, Qoal^'^ 53.14 meV, 

0.067me 



S[ ^ai = 0.4916, m* 



and K = 12.4. 



The anisotropic response of the electron system to the 
two different polarizations of the photon field is caused 
by the confinement energy in the ^/-direction, HQq = 1.0 
meV, beeing larger than both the photon energy, hu = 
0.4 meV, and the lowest excitation in the x-direction, 
AEx ~ 0.2 meV. At the low energy we are observing 
the properties of the system here it is much "harder" 
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